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c3 ! Abstract 



A phenomenological theory of fermion masses and mixings is constructed 

^ ' within the framework of a four- family symmetry. It is found that the 

^ ; most favored set of relevant CKM elements are |T4s| ~ 0.222, \Vcb\ ~ 0.044, 

O ■ \yuh/Vd,\ ~ 0.082, \Vud\ ~ 0.974, |F,,| ^ 0.9736, \V^d\ ~ 0.224 with Bk ^ 0.8. 

0^ ' The top quark mass is predicted to be 258 GeV at 1 GeV with its physical 

Q^' mass approximately equal to 153 GeV. The Majorana scale associated with 



the fourth neutrino is bound from above to be 6.4 TeV. 
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Despite the impressive agreement of the standard model with experiment, it is clear that 
questions such as the origin of fermion masses and mixing angles cannot be answered solely 
within the framework of the standard model. At the present time, one does not know at 
what energy scale ( called Family scale here) lies the solution to the mass problem. There 
is a some belief that such a scale might be very near the Planck mass where it is hoped 
that all interactions (including perhaps gravity) are unified. This might be the case. On the 
other hand, it is also possible that the Family scale might not be too much higher than the 
electroweak scale. We shall explore this possibility below. 

It is a fact that not only do we have inter-generation mass splitting but we also have 
intra-generation mass splitting. This last splitting breaks explicitely the so-called " custodial" 
SU{2) symmetry (in the simplest version of the standard model) which guarantees the p 
parameter to be equal to unity at tree level. What is the relationship between the " custodial" 
breaking and the behaviour of the fermion mass matrices? One first notices that, if the mass 
matrices for the Up and Down sectors were identical (in form and magnitude), the CKM 
matrix []I| would become the unit matrix and all mixing angles would vanish. The fact that 
the CKM mixing angles are non-vanishing and small is an indication that the fermionic 
"custodial" breaking should be non- vanishing and small as well. By "small" we mean that 
the contribution to the p-parameter is small. We shall use this hint in the construction of 
our model. 

To address the mass problem, it is obvious that one has to go beyond the minimal (three 
generations) standard model. We propose in this note a model with four generations and 
a family symmetry p|. This (broken) family symmetry gives rise to mixing terms in the 
fermion mass matrices at one-loop level. One of the motivations for having a fourth family 
here is the fact that, in our model, the main contribution to the off-diagonal elements of 
the mass matrices come from one- loop diagrams with a scalar exchange, and large Yukawa 
couplings are preferable to guarantee that these elements are non-negligible. The nature of 
fermion masses and mixing in our framework is linked to the nature of the scalar vacuum 
expectation values and resulting scalar mixing. We shall show that the mixing terms in the 



quark and lepton mass matrices are intrinsically related. As a consequence, we shall also 
show how a rare decay process such as Kl -^ /xe can set an upper bound on the Majorana 
mass scale. 

To set the stage for the construction of the mass matrices, we now list our assumptions. 
They are: 

1) There is a vector-like S'[/(4) family (gauge) symmetry. 

2) The combined generation-electroweak symmetry is S'f/(4) ® SU{2)l ® SU{2)r ® 
U{1)b~l- The existence of SU{2)ji in our context is linked to the assumption that the 
family symmetry is vector-like. 

3) There is a Higgs field, $ = Tj$* with i = 1, . . . , 15, which transforms as (15, 2, 2, 0) and 
another one, 0, which transforms as (1, 2, 2, 0) under SU{4:) (g) SU{2)l ® SU{2)ii ® U{1)b-l- 
(We are ignoring here scalars which are electroweak singlets.) 

4) The left and right-handed quarks transform respectively as (4, 2, 1, 1/3) and 
(4,1,2,1/3) while the left and right-handed leptons transform as (4,2,1,-1) and 
(4, 1, 2, -1). The charged generator is Q = T^l + ^3^ + {B - L)/2. 

5) There exists a Yukawa coupling of the form 

Cy = qiihi^ + /i2$ + Kcj) + /i20)gi? + h.c. 

+lL{h^ + h^ + l[(j) + 124>)Ir + h.c. (1) 

where the notations used for quark and lepton fields are self-explanatory and where ($, 0) = 
(Tir2<l>^V2,r20V2). 

6) The Higgs fields have vacuum expectation values (VEV) 

($) = (($3)/v^)t/m(7(l, -1, 0, 0) + i{^s)/V6)dtagil, 1, -2, 0) + {{^i,)/Vl2)dtagil, 1, 1, -3), 

(2a) 

{{p) = diag{K, k), (2b) 

where 



($3,8,15) = diag {V{3,8,15} , ^3,8,15})> (3) 

The various VEV's will in general be complex and can be parametrized generically as ve''^. 
7) The breaking of the family symmetry leads to the following mass matrices for the 
scalars which connect the fourth family to the other three and to itself, and which are 
•^1,4 = ($9 T i$io)/V2, <^2,4 = ($11 T i$i2)/v^, $3,4 = ($13 T ^$14)7^2, and $44 = $15, 



M' 



' 1 a 0^ 



a 1 b 
6 1c 
c 1 



(4) 



for the neutral sector, $°, and 



1 a' ^ 



M' 



(5) 



/ 



a' 1 b' 

b' 1 c' 

c' 1 

for the charged sector, $'^. The diagonal elements are (14), (24), (34) and (44) from top to 
bottom. The six coefficients a, . . . ,c' will determine the amount of mixing among different 
generations. It is found that the primed and unprimed coefficients are not too different from 
each other, which means that the amount of scalar "isospin" breaking is small. We shall 
assume that the mixing among other generation-changing scalars, if it exists, is negligible 
compared with the above and hence the only dominant off-diagonal elements of the mass 
matrices are those generated by the above scalars. 

8) Since $ transform as (2, 2) under SU{2)l ® S'f/(2)/j, each component $* takes the 
form 



$« 



^ 0?^ 02^^^ 



V 



-0r^ 0^ 



(6) 



In the computation of the one-loop contribution to various off-diagonal elements, both 0| 
and 02 contribute. We will assume that they both have mass mixing of the form of Eqs.(4,5). 



Furthermore, we will assume that the coefficients a, . . . , c' of the two sectors are proportional 
to each other, with a common proportionality coefficient denoted by y. 

Flavor-changing neutral currents can be suppressed provided the masses of the family 
gauge bosons and mixing parameters in both the left-handed and right-handed quark "ro- 
tation" matrices are appropriately chosen. Similar considerations apply to the scalar sector. 
It is beyond the scope of this paper to present such an analysis. We shall assume it can be 
done and shall concentrate here on the construction of the mass matrices. 

Various mixing parameters and phases in the fermion mass matrices are explicitely com- 
puted in terms of parameters of the basic family symmetry. Once the mixing parameters 
are fixed, the phases are completely determined. These same parameters enter in the mass 
matrices of both quark and lepton sectors so that the mixing angles and phases of the two 
sectors are found to be related. 

The mass mixing among the scalars are assumed to be such that only the following 
mixings occur: ^44 — ^43, ^43 — 042, 042 — 04i giving rise to a particular type of symmetric 
mass matrices. In a way, fermion mixing can be seen as a direct consequence of scalar 
mixing. If we recall Eq.(2), the breaking of the family symmetry is assumed to be such 
that all of the family changing scalars have vanishing vacuum expectation values (a familiar 
assumption). In this case, there is no direct fermion mixing at tree level coming from the 
symmetry breaking. On the other hand, a general gauge- invariant scalar potential will 
yield, after symmetry breaking, various mass mixings among the scalars. Couplings of the 
fermions to these same scalars will give rise to fermion mixing beyond the tree level in the 
mass matrices. The Yukawa couplings that enter these radiative corrections are much larger 
than the gauge couplings and, because of that feature, we shall ignore the contributions 
coming from the family gauge bosons. 

The diagonalization of Eqs.(4,5) can be done analytically in a straightforward manner. 
Each of these four scalar mass eigenstates can now couple the fourth generation to all the 
other three, with the strength given by the elements of the eigenvectors. The trivial details 
are given elsewhere. The only remark one would like to make here is the fact that the mass 



scale M^ mentioned above do not enter the calculations given below because the results 
there (inside the logarithms) are given in terms of mass ratios. 

At the one-loop level, the off-diagonal elements of the fermion mass matrices are gen- 
erated by the exchange of the above four scalars and are finite In our model the twenty 
four off-diagonal elements of the up- and down-symmetric mass matrices of the quark and 
lepton sectors can be computed in terms of the six parameters (a, . . . , c'), the four fourth- 
generation masses, and a phase difference (to be specified below). Let us first concentrate 
on the following elements in a typical mass matrix: 4-3, 3-2, 2-1. For the up sector (for 
quarks), the one-loop contributions are given by 

Mg = hi e''' |M° I {\"xe-'^ - -)Ll (7a) 



Mg = hi e''^ |M° I {Xxe-'^ - -)Lo, (7b) 



Mg = hi e'^' |M° I {X'xe-'^ - -)L[„ (7c) 

where r = /i2/^i and x = \M^\/\Mjj\, the ratio of the fourth-generation quark masses. For 
the down sector, one has 

Mg = hi e'^' |M° I (A"e*^ - rx)Ll (8a) 

Mg = hi e'^'' |M° I (Ae^^ - ra:)Lo, (8b) 

Mg, = hi e'^' |M° I (AV^ - rx)L'Q. (8c) 

Here, Lq = L2 — L, . Assumption (8) gives L2 = r'^L, (here y = r^) so that 
Lq' = (r^ — l)Li . The Lq's are the one-loop contribution to the off-diagonal elements 
coming from neutral scalars. The parameters A*^''"^ represent the contribution of the charged 
scalars relative to the neutral ones. They are defined for the case when r 7^ 1. Assumption 
(8) gives Lp''"^ = r^Lt^''"\ One then has A^'") = {r/{r^ - 1))(L+^''"V^?^''"^ Notice that 



the various L's can be explicitely computed using the mass eigenstates and eigenvectors of 
Eqs.(4,5). They are typically of the form: {Aln{mi/m2) + Bln{m^/ m^) / IQn"^ , where A and 
B are some combinations of the elements of the eigenvectors. The parameter A is defined 
as A = 5i — ^2 where m^ = 6**^^1772^1 and m^, = e^^'^\m^j^\. 81^2 are functions of the phases of 
the complex VEV's. 

The parameters A's and Lq's depend entirely on the scalar sector and that is the reason 
why they appear in both Eq.(7) and Eq.(8). This will also be the reason why they also appear 
in the lepton mass matrices. There one just has to make the replacements: hi^2 -^ h,2, 
r ^ ri = h/h and x — ;► X; = |^£;|/|^Ar| in Eqs. (7,8). 

Notice that we seem to ignore elements such as M42 and M41 which can be computed 
at the one-loop level and other off-diagonal elements which arise at higher order. The 
reason is simply that they are numerically small compared with the above elements. Within 
the precision considered in this paper, they are found to be not so important. A more 
detailed and precise analysis to be carried out in a subsequent work will include these extra 
corrections. It is however important to state here the fact that these extra terms are perfectly 
calculable in terms of known parameters of the model. 

Let us now examine Eqs. (7,8) in the case when the "isospin" breaking parameter 
r = h2/hi becomes unity, i.e. no "isospin" breaking. A look at the diagonal masses obtained 
from Eqs. (1,2,3) reveals that when r = 1 or equivalently h2 = hi, one obtains rriij = m'jj, 
m^ = m^, 777° = 777°, and 777^ = m^. Now 777^ = 777^ means that the parameter x = 1 and that 
Si = 62 giving A = 0. Eqs. (8,9) then tell us that M^^^ = Mf^, M^^ = M^^, and M^^ = M^^. 
(One can make similar statements for the other off-diagonal elements.) This means the the 
mass matrices for the up and down sectors are identical and so are the matrices Uu and Uu 
which diagonalize them. This now means that the generalized CKM matrix Vckm = U^^Ud 
is equal to the unit matrix. All mixing angles vanish. This is precisely the point that we have 
mentioned above. This feature is true regardless of how much "isospin" breaking there is 
in the scalar sector, i.e. regardless of how much mass splitting there is between the charged 
and neutral scalars. In our model it turns out that the magnitude of the CKM elements is 
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a result of the interplay between the "isospin" breaking term of the Yukawa sector, r, and 
that of the scalar sector as defined by the parameters A'-''''^ 

In order to carry out a numerical study of mass eigenvalues and CKM matrix elements, 
it is found to be more convenient to parametrize the off-diagonal elements of the mass 
matrices in a slightly different way although in principle one can just use directly Eqs. 
(7,8). One just has to vary the parameters described above in order to fit the generalized 
CKM matrix elements and the quark masses. This is a perfectly well-defined task albeit 
a very time-consuming one. To narrow down the range of values, we are guided by the 
hierarchical nature of the masses and the sizes of the CKM matrix elements. We shall 
parametrize the various M's as: M.^^ = cuXfjC'^^u^ -^32 = f^t^u^^^'^ , -^21 — rtXfje'^^u ^ for 
the up sector, and Mf^ = cz)A|,e*'^D M§2 = rbXce'^^^, M.21 = rbXj^e'-^'D , for the down 
sector. The previous quantities are the elements of A4u and A4d where Mu = \M^\e'^^'^A4u 
and Md = \M^\e'^^'^A4d■ The parameters rt and r^ are defined as r^ = |m^|/|mf^| and 
Tf, = \m^\/\mfj\. The above phases are defined as tan5jj = XxsinA/{l/r — XxcosA) and 
tanSo = XsinA/{XcosA — rx) and similarly for the primed and double-primed quantities. 
The phases ^1^2 can be absorbed in a redefinition of the up and down quark fields. The 
parameters AJy'J are defined in terms of those of the fundamental theory via Eqs. (7,8). In 
total the ten parameters A[/, Ad, cjji cd, Su'jj are computed interms of seven parameters 
r,X^''"\ Lq/Lq, Lq/Lq, and A. To complete the picture, one has to specify the diagonal 
elements coming from the other three generations for both up and down quarks. In general 
they have arbitrary phases and magnitudes. Although the phases can be similar to those of 
the fourth generation, we shall treat them as free parameters. It turns out that the attractive 
possibility of having all phases approximately equal to each other can actually be realized 
in our model. 

After factoring out |M^|e*'^^'2 in the up and down mass matrices, the first three diagonal 
elements can be written as r„e*"% TcC*"^, r^e*"^ for the up sector and r^^e*'^^, r^e*^^, rf,e*^^ for 
the down sector, where ai = 5„ — (5i, ^2 = ^c — ^i, 0:3 = 6t — Si, and Pi = 5^ — ^2, P2 = ^3 — ^2, 
Ps = Sb — 62- Here r^ = |m^|/|m^|. Aiu and Aio which are symmetric matrices can be 



made real by an appropriate redefinition of the quark phases provided ai = 2{S'{j + S[j — Su), 
a2 = -2{6'/j - 6u), as = 26'{j, A = 2(5^ + 5^ - 6d), ^2 = -2(5^ - 6d), /^s = 25^. Assuming 
that Su,Sc, etc... satisfy the previous "reahty" condidtions, it will be seen that the only 
CP phases which enter the generalized CKM matrix are A = 6u — So, ^' = S'u — 6'jj, and 
A" = S'Ij — 6'1). Phenomenologically, Si/,6d, etc... are found to be either close to (but not 
equal to) 0° or 180° so that Su,Sd, etc... are almost equal to 61 and 62- Such a possibility is 
rather attractive: all quark masses have roughly similar phases. 

To bring Aiu and TWd to a real form, we redefine the left-handed quark phases by a 
diagonal matrix of the form Qu = diag{l,e~'^'^^ ,e~^'^2 , e~*'^3 ) for the up sector and Qo = 
diag{l,e~^'^i ,e~^'f'^ , e"*'^^ ) for the down sector, where 0i = 25" + 6' — 25, 02 = S' — 5, 
03 = 6" + 6' — 6 with the appropriate subscripts for up and down. (There is also a redefinition 
of the right-handed quark phases but it is irrelevant for Vckm and we shall not discuss it 
here.) The eigenvalues of the real, symmetric 4x4 matrices denoted by Ai§ and Ai^ can 
be computed numerically. These eigenvalues are in turn used to construct the eigenvectors 
which are then used to obtain the matrices Ru and Rd which diagonalize Ai^ and Ai^ 
respectively. The generalized CKM matrix is then defined as Vckm = Qu^Rjj^RdQd- A 
full account of our analysis will be given in a separate publication. We shall illustrate 
our model with a few examples in this letter. A few typical elements are Vus = {^%/\''^s\ — 

(f„/A3,)e^(2A"+A'-2A) +0(^3)1/^^^^ V^^ = (~J^2^)|;^-lg.(A'-A) + (;^^/^^^^)g.(2A"+A'-2A) + 

0(A2)}/iV„5, Kb = {-(ir,/Az,)e^(2^"+^'-2^) - (A^/ir,)e^(^'-^) + 0{X')}/N,t,. All other 
elements of the generalized CKM matrix are easily computable and they will be given in 
an extensive version of this manuscript. The various A^'s are normalization factors. The 
f's are related to the ratios of the eigenvalues and r^ or r^. The quantities i^j's are Ki = 
1— fj — (c(7,£)Af^£))^/(r(_fe(l — Tj)) with i = u,c,b. One can immediately notice that the mixing 
of the "light" generations with the fourth generation creeps into Vuh and Vet via various Kj's 
which contain Cu or C£, which enter ^^43. This is because of the way the fourth generation 
mixes with the third one in our mass matrices. The CKM elements which do not involve 
the third generation-we are mainly concerned with the 3x3 sub-matrix here-are relatively 



insensitive to the presence of the fourth generation. In this sense, the physics of the third 
generation (in particular B-physics) indirectly probes the existence or non-existence of a 
fourth generation. Recall that the fourth generation indirectly manifests itself through the 
parameter x which enters in all elements of the mass matrices. It is however the resulting 
(dominant) mass mixing with the third generation that manifests itself most visibly in the 
CKM elements involving the third generation. 

An extensive numerical analysis is underway. We shall give here some preliminary re- 
sults. The inputs are given in the form of M(ri,r2,r3, 1) for the masses. We have |^ 
\M^\ = ITeV, X = \M%\/\M^\ = 0.98, lTeV{-.8 x 10-^ .00213, .54) for the up sec- 
tor and 0.98re1/(.596 x 10-^ -0.952 x 10~^ .0062) for the down sector, Xu = 0.064, 
Xd = 0.15, C{/ = 112.5, cn = 0.325, r = 1.4637. The previous values not only deter- 
mine the mass eigenvalues but they also fix the phases in the following way. By equating 
the two ways of writing M43, M32, and M21 for the up and down sectors and by tak- 
ing the absolute value of these elements, we can fix the values of A' and A" once A and 
A are given. Let us recall that the phases Su, 5o and their primed and double-primed 
counterparts are given in terms of x, r and A, A', and A" respectively. Since the CKM 
elements depend on the phase differences A, A' and A" defined above, we have the fol- 
lowing four sets of values: 1)(1.47726, 1.262553, 1.4592) giving (7.962°, -181.124°, -1.299°), 
2)(1.477385, 1.26252, 1.459415) giving (6.508°, -180.92°, -1.0631°), 

3)(1.4775, 1.2625, 1.45957) giving (5.0566°, -180.715°, -0.8269°), 

4)(1.47753, 1.262495, 1.45965) giving (4.334°, -180.613°, -0.7088°), where in (1) to (4) the 
first and second sets correspond to (A, A', A") and (A, A', A") respectively. The results are:l) 
IK.I = 0.22, \V4 = 0.0476, iKb/K^I = 0.0834, |K,| = 0.9739, \VJ = 0.222, iKdl = 0.9745, 
Bk = 0.555, 2)|Ks| = 0.221, |Kfe| = 0.0459, |WKfe| = 0.0827, |K.| = 0.9737, \VJ = 
0.2228, iKdl = 0.9743, Bk = 0.6385, 3)|K.| = 0.222, |Kb| = 0.0444, |Kfe/Kfe| = 0.082, 
\Vcs\ = 0.9736, \Vcd\ = 0.2237, |Kd| = 0.9741, Bk = 0.794, 4)|Ks| = 0.2225, \Vcb\ = 0.04378, 
\Vub/Vcb\ = 0.0816, \Vcs\ = 0.9736, |Kd| = 0.224, |Kd| = 0.974, Bk = 0.913, for the four 
cases. In all four cases, the mass eigenvalues are (at IGeV) mjj = 1.28TeV^, nit = 258G'eV^, 
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mc = l.SGel/, rriu = 5.16MeV, tud = 980Gel^, m^ = 6.15GeV, m, = U7MeV, 
rrid = 8.6MeV. The values of B^ were obtained using the experimental value of e namely 
0.00227. In this preliminary analysis the physical top quark mass is mf ^ 153GeV. A 
better way to present the above results would be to use a plot of the CKM elements 
versus B^ but the trend can already be seen from the above numbers. The latest lat- 



tice computation of the Isgur-Wise function 0] give {VcblyTB/lASps = 0.038^2^3 [H 



\Vcb\^TB/l.53ps = 0.044 ± 0.005 ± 0.007 §, and Bk = 0.825 ± 0.027 ± 0.023 using stag- 
gered fermions 0. Taken at face value, Bk tends to favor \Vus\ ~ 0.222, \Vcb\ ~ 0.044, 
\Vub/Vcb\ ~ 0.082, \Vud\ ~ 0.974, |K.| ~ 0.9736, \Vcd\ ^ 0.224. Further experimental and 
theoretical efforts are clearly needed to pinpoint these values. 

One might also ask about the sensitivity of the top quark mass prediction to the value of 
the fourth-generation mass. Our preliminary results indicate that, for the fourth-generation 
quark mass not too far from 1 TeV, the result is not very sensitive to the precise value of 
that mass since having fixed Xu and Xn, one has rrit ^ rric/Xfj and therefore rrit is more 
sensitive to Xu and rric than mij. A more detailed investigation of the dependence on rriu 
for a larger range of values will be carried out separately. The results presented here are for 
mu around 1 TeV. 

The lepton sector is identical in form to the quark sector and the results can be obtained 
with the replacement hi 2 -^ /i,2- A rough estimate gives A^; ~ lAXD{m^/m^){m^j^/mfj)^ 
and Xn ~ 1AXe{Xu/Xd)- Here m^^ is the Dirac bare mass of the fourth-generation neutrino. 
Since \Vus\ ^ Xu + Xd, the leptonic version is iKe^l ^ Xn + Xe- A back of the envelope 
estimate gives BR{Kl -^ fie) ^ {X^ + XeYBR{Kl -^ fifi). With BR{Kl -^ fifi) = 
7.3 X 10-s and BR{Kl ^ /le) < 9 x 10-i\ we find m% < 451 GeV. For the sake of 
estimate, the mass eigenvalue is assumed to be similar to that obtained for mu, namely a 
fator of 1.2 of its input, giving rriNphys < 541 GeV. Assuming a see-saw mechanism of the 
form m^^ = mNphysl -^ ^-^^ requiring that m,^^ > 46 GeV, one obtains M. < 6.4 TeV. 
A more detailed calculation would presumably gives a bound on the Majorana scale (at 
least for the fourth generation) not too far from the previous value. Also it is not hard to 
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arrange the masses of the other three neutrinos to be much hghter than the fourth one. A 
host of interesting phenomena might be studied such as an intriguing possibihty of direct 
CP- violation in r-decay, presumably in some kind of r factory. 

This work was supported in part by the U. S. Department of Energy under Grant No. 
DE-A505-89ER40518. 
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